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Abstract We compare two approaches to nonequilibrium thermodynam- 
ics, the two-generator bracket formulation of time-evolution equations for 
averages and the macroscopic fluctuation theory, for an isothermal driven 
diffusive system under steady state conditions. The fluctuation dissipation 
relations of both approaches play an important role for a detailed compari- 
son. The nonequilibrium Helmholtz free energies introduced in these two ap- 
proaches differ as a result of boundary conditions. A Fokker-Planck equation 
derived by projection operator techniques properly reproduces long range 
fluctuations in nonequilibrium steady states and offers the most promising 
possibility to describe the physically relevant fluctuations around macro- 
scopic averages for time-dependent nonequilibrium systems. 
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1 Introduction 

There appears to be a widespread believe that time-dependent far-from- 
equilibrium systems are too complicated to be dealt with and that one should 
hence develop a thermodynamic framework for nonequilibrium steady states 
first. With a strong focus on phenomenological thermodynamics, we just men- 
tion the fundamental attempts to establish a general theory of nonequilibrium 
steady states made by Oono and Paniconi (1998) [41], Sasa and Tasaki (2006) 
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[62], Bertini et al. (2002-2009) [6 l l7 l [8 l l9] . and Taniguchi and Cohen (2007) 
[64] . This hst of efforts is by no means complete, but it gives a flavor of the 
variety of different approaches; many further references can be found in the 
above papers. Most of these groups look for support for their phenomeno- 
logical thermodynamic approaches from statistical mechanics. Among the 
important insights into the statistical mechanics of nonequilibrium steady 
states, we mention the work of Derrida et al. [TSllTB] on the fluctuations in 
an exactly solvable model of a driven diffusive system, fluctuation theorems 
[2T|l23ti22j . and the fundamental concept of SRB measures (Sinai, Ruelle, 
Bowen [63 ] ll0 1 fTT 1 [60]) defined on the attractors of chaotic systems (see also 
the reviews [TSlleTl leS]). 

While the relative simplicity of steady state systems seems to be indis- 
putable, there may also be some reasons to believe that looking in a suffi- 
ciently abstract way at the general problem of nonequilibrium time evolution 
might help to clarify the structure of nonequilibrium thermodynamics. First 
of all, the formulation of time evolution can be given in terms of illuminat- 
ing geometric structures for generating trajectories. Second, in the evolution 
equations one can easily separate reversible and irreversible contributions, 
where the hallmark of reversibility is the possibility of a Hamiltonian for- 
mulation. Third, the fundamental thermodynamic concepts of energy and 
entropy reappear naturally in the leading parts of nonequilibrium dynamics 
as generators of reversible and irreversible evolution, respectively. Fourth, in 
time-dependent local field theories one can easily separate bulk from bound- 
ary effects, which is difficult, if not impossible, in a steady state system be- 
cause the influence of the boundary conditions penetrates through the entire 
system. Fifth, and perhaps most importantly, we have the machinery of the 
projection operator formalism [66 p 38 1 f37 p 59 p 25j . which relies on a clear sepa- 
ration of time scales. Separation of time scales is a trump card gamed away in 
the theory of steady state systems. Detailed elaborations on all these issues 
can be found in the textbook which summarizes a stream of developments 
that started with three pioneering letters [29 1 139 1 126] in 1984, was concretized 
into the clarifying textbook |4], and formulated as a general framework in 
the papers [27,.,.55,44j. 

This paper is organized as follows. We flrst summarize the two-generator 
bracket formulation of the general evolution equations for nonequilibrium 
systems (Sec. [5]) and we then illustrate the ideas of this framework by deriv- 
ing the equations for a driven diffusive system (Sec. [3]). In SeclH macroscopic 
fluctuation theory is presented in the context of the driven diffusive system. 
The different approaches to nonequilibrium thermodynamics are compared 
in Sec.ini which leads to the distinction between thermodynamics of the first 
and second kind based on averages and fluctuations, respectively. In Sec. [6l 
fluctuation dissipation relations are used to establish deeper connections be- 
tween the two-generator bracket formulation and the macroscopic fluctuation 
theory. We give an explicit example of a calculation of long range correlations 
in a nonequilibrium steady state within the phcnomcnological two-generator 
bracket approach. A summary and discussion conclude the paper (Sec.[7|). 
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2 GENERIC 

Time-evolution equations for nonequilibrium systems have a well-defined 
structure in which reversible and irreversible contributions can be specified 
separately. In particular, the reversible contribution is generally assumed 
to be of the Hamiltonian form and hence requires an underlying geomet- 
ric structure (a Poisson bracket) which reflects the idea that the reversible 
time evolution should be "under mechanistic control." The remaining irre- 
versible contribution is generated by the nonequilibrium entropy by means 
of a dissipative bracket. The nonequilibrium energy and entropy landscapes 
are introduced through their roles as generators of reversible and irreversible 
dynamics in the space of nonequilibrium variables; they are associated with 
the evolution of averages. It has been shown in Sec. II. B. 5 of [27] how equi- 
librium thermodynamics in its familiar form arises from these generators of 
time evolution. 

Our discussion is based on the general equation for the nonequilib- 
rium reversible-irreversible coupling (GENERIC) for the time-evolution of 
nonequilibrium systems p7[f55ll46| . If A is an arbitrary observable, that is, 
a sufficiently regular real-valued function or functional of a set of variables 
X required for a complete description of a given nonequilibrium system, the 
time evolution of A is given by 

±- = {A,E} + [A,S]. (1) 

The observables E and S generating time evolution are the total energy and 
entropy, and {•, •} and [•, •] are Poisson and dissipative brackets, respectively. 
The bracket of two observables A and B is another observable with a linear 
dependence on A and B (a more complete characterization of Poisson and 
dissipative brackets is given in Eqs. ([I])-® below). The two contributions 
to the time evolution of A generated by the total energy E and the entropy 
S in Eq. ([T]) are referred to as the reversible and irreversible contributions, 
respectively. Equation ([1]) is supplemented by the complementary degeneracy 
requirements 

{S,A} = 0, (2) 

and 

[E,A]=0, (3) 

which hold for all observables A. The requirement that the entropy is a de- 
generate functional of the Poisson bracket expresses the reversible nature of 
the first contribution to the dynamics: the functional form of the entropy is 
such that it cannot be affected by the Poisson bracket contribution to the 
dynamics, no matter which observable A is used as a generator of reversible 
dynamics. The existence of degenerate observables is a hallmark of coarse 
graining because the Poisson bracket associated with the symplectic struc- 
ture of atomistic equations is non-degenerate and hence does not allow for the 
existence of an entropy on the purely reversible atomistic level. The require- 
ment that the energy is a degenerate functional of the dissipative bracket 
expresses the conservation of the total energy by the dissipative contribution 
to the dynamics in a closed system. 
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For completeness, we give the defining properties of Poisson and dissipa- 
tive brackets. The Poisson bracket possesses the antisymmetry property 

{A,B} = -{B,A}, (4) 

and satisfies the product or Leibniz rule 

{AB,C} = A{B,C] + B{A,C], (5) 

as well as the Jacobi identity 



{A, {B, C}} + {B, {C, A]] + {C, {A, B}} = 0, (6) 

where A, _B, and C are arbitrary observables. These properties are well-known 
from the Poisson brackets of classical mechanics, and they express the essence 
of reversible dynamics. The Jacobi identity ([6]), which is a highly restrictive 
condition for formulating proper reversible dynamics, expresses the invari- 
ance of Poisson brackets in the course of time (time-structure invariance) . 

The dissipative bracket satisfies the symmetry condition (for a more so- 
phisticated discussion of the Onsager-Casimir symmetry properties of the 
dissipative bracket, see Sections 3.2.1 and 7.2.4 of [46]) 

[A,B] = [B,Ai (7) 

and the non-negativeness condition 

[A,A]>Q. (8) 

This non-negativeness condition, together with the degeneracy requirement 
(121), guarantees that the entropy is a nondecreasing function of time, 

f = [^,^]>0. (9) 

The condition ([5]) may hence be regarded as a strong formulation of the 
second law of nonequilibrium thermodynamics. 

In practical calculations, it is often convenient to formulate GENERIC 
in terms of Poisson and friction operators instead of brackets [551l46j . This 
is the analogue of formulating Hamilton's equations of motion in terms of a 
symplectic matrix rather than in terms of Poisson brackets. More precisely, 
one generally writes 

5A5B 



and 



where L is the Poisson operator and M is the friction operator. In case that 
L and M are differential operators, boundary terms need to be discussed 
separately, as illustrated in the subsequent section on diffusive systems. The 
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Leibniz rule ([5|) then follows automatically. The time-evolution equations for 
the system variables x implied by Eq. ([T]) can be expressed in the form 

dx ^6E (12) 

dt 6x Sx 

If there is a constant Tq with dimensions of temperature, then we can intro- 
duce the Helmholtz free energy, 

F = E- TqS, (13) 

and rewrite Eq. ([T^) in terms of the single generator F, 

dx SF 1 ^SF 

where the mutual degeneracy requirements ^ and ^ have been used. Dou- 
ble generator {E, S) and single generator (F) formulations of nonequilibrium 
thermodynamics have been compared in great detail [T9 j l20 l l3]. 

The safest case of a constant temperature arises in quantum gravity where 
To is the Planck temperature, that is, a constant of nature [52]. In the present 
context, however, the idea is to consider isothermal systems. Although very 
common, the assumption of isothermal conditions is subtle and may even 
obscure fundamental discussions of nonequilibrium systems. As entropy pro- 
duction is a key feature of nonequilibrium systems, it would be preferable to 
include heat conduction into the discussion. Steady states are naturally as- 
sociated with nonuniform temperature profiles and the transport of entropy 
through the boundaries of an open system. The treatment of the required 
boundary conditions within the GENERIC framework, which is crucial for 
the discussion of driven systems, has been discussed in a number of recent 
publications [171I5I53..54J. In the context of SRB measures, isothermal con- 
ditions are achieved by including a thermostat into the equations of motion, 
which might be considered as undesirable for an attempt to clarify the con- 
ceptual foundations of nonequilibrium thermodynamics. 

The structure of GENERIC can be obtained by the projection operator 
method [44 l l45 l l46] . This technique relies on nonequilibrium statistical en- 
sembles, in particular, generalized microcanonical and canonical ensembles, 
and projectors on the spaces of slow and fast variables. As a result, one 
obtains practical recipes for calculating the GENERIC building blocks by 
means of statistical mechanics [49]. In this paper, however, we focus on the 
phenomenological approach to nonequilibrium thermodynamics. 



3 Driven Diffusive System 

To illustrate the GENERIC framework, we derive the evolution equation for 
a diffusive system. As our only system variable x, we choose the mass den- 
sity field X = p{r) of the diffusing species on a domain A. We further assume 
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isothermal conditions at temperature Tq. For this assumption to be mean- 
ingful, we consider an athcrmal hard-sphere system for which the internal 
energy density is entirely of kinetic origin and hence of the ideal gas form, 

3 /cbTo , , 

where fee is Boltzmann's constant and m is the mass of a spherical particle. 

An important observation is that the derivative of e(p) with respect to p is 
constant. The entropy density of a hard-sphere system is given by 



s{p) = —pin 

m 



rr3/2 Ro{p) 

:; 



(16) 



where Ro{p) is a given function of p. For an ideal gas, that is, in the limit of 
vanishing particle radius, Ro{p) is a constant with proper dimensions (to be 
constructed with the help of Planck's constant). The derivative of s{p) with 
respect to p at constant Tq is associated with the chemical potential p per 
unit mass according to 

p 1 df 1 de ds ds 3 /cb 

Th^ Todp^ Thdp~ dp^ ~dp^ 2^' ^ ' 

where the Helmholtz free energy density for our isothermal system is given 

by 

f{p) = e{p)-Tos{p). (18) 

The total energy and entropy are obtained as the integrals of their den- 
sities over the domain A, 

E= f e(p{r))cPr, S = [ s{p{r))d^r, (19) 

J A J A 

and, for these simple functional, the functional derivatives of E and S with 
respect to p{r) are simply given by the partial derivatives of the respective 
densities. 

As a next step, we need to introduce the Poisson and friction operators 
L{x) and M{x). For hydrodynamic systems, the only reversible effect is con- 
vection. As we consider a purely diffusive system and we have eliminated 
the velocity field from the description, the mass flux is treated as entirely 
diffusive, that is, L{x) = 0. For the friction operator, we choose the local 
diffusion operator 

^ = -|:-^°^-|:' (2°) 

which, for a positive semidefinite symmetric tensor S, is also positive scmidef- 
inite (as can be shown after an integration by parts and after ignoring bound- 
ary terms for the purpose of deriving bulk evolution equations). These prop- 
erties are even more obvious in the bracket notation 



dr Sp J \dr 6p 
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The degeneracy requirement ([3]) for this dissipative bracket is satisfied be- 
cause 5E/5p is constant for our hard-sphere system. We have now specified all 
thermodynamic building blocks of the GENERIC framework and are ready 
to write out the evolution of our diffusive system. From the fundamental 
equation (|12p we obtain the evolution equation 

^ = A . s • ^ (22) 

dt dr dr 

The gradient of the chemical potential is the natural driving force for diffusive 
mass transport in an isothermal system. With the help of the explicit form 
(|16p of the entropy the driving force can alternatively be formulated in terms 
of the density gradient, 

dp ^ d_ kBTp R2{p) - . ^ ^ _^ . N ^ /23-v 
dt dr m p dr dr dr' 

where the function i?2(p) results from first and second order derivatives of 

p2 dRoipY 



R2{p) = l-j- 
dp 



(24) 



_Ro{p) dp 

and the second part of Eq. (I^H)) is merely a definition of the diffusion tensor 



D{p)^^^S. (25) 

m p 

Even for a constant tensor S, a complicated dependence of the diffusion 
tensor on p arises. For the ideal gas, we find i?2 = 1 and D (x 1/p. 

The rate of change of the total energy in the system is given by the chain 
rule gZ], 

dE SE^ JS ^ SfcBTo f ^ dp 2 



dt 5x Sx ^ ' ^ 2 TO Jgy^ dr ' ^ ^ 

where n is the normal vector on the boundary dA of the domain A pointing 
from the inside to the outside of the system. The boundary term is a result 
of the integration by parts required to go from Eq. (^0]) to Eq. (PT|) . As 
the dissipative bracket [E, S] vanishes according to Eq. ^ , the energy of 
our isothermal hard sphere system changes according to the mass flux into 
the system. For steady state systems, the total mass flux must be zero. We 
are interested in the nonequilibrium steady state systems arising for a given 
nonuniform chemical potential on dA. 

With the bracket in Eq. ((2T|) . wc have introduced a dissipative mass 
flux. Such a possibility has been debated controversially in the context of 
the full Navier-Stokes-Fourier equations of hydrodynamics. In a classical pa- 
per, Dzyaloshinskii and Volovick |T7] proposed the inclusion of a dissipative 
mass flux into the hydrodynamic equations. Starting from a modified kinetic 
theory, Klimontovich [3D] arrived at the same suggestion. More recently, a 
dissipative contribution to the mass flux was re-introduced and forcefully 
promoted by Brenner [12J, whose work stimulated significant interest and 
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controversy in the physics and fluid dynamics communities. Brenner's work 
motivated thorough investigations on the thermodynamic admissibihty of a 
dissipative contribution to the mass flux '46',f2l[Tl l48l[50] , which demonstrated 
that the idea fits naturally into the GENERIC framework and into standard 
linear irreversible thermodynamics. However, these investigations focused en- 
tirely on nonequilibrium thermodynamics and neglected other, equally im- 
portant considerations, such as the local conservation of angular momentum. 
Such additional criteria were considered earlier in ^31] and corrobated the 
original (not rigorously justified) statement of Landau and Lifshitz that a 
dissipative contribution to the mass flux cannot exist (see footnote at the 
end of Sec. 49 of [35]). A concise summary of the current state of the dis- 
cussion can be found in the comment [35] to the letter [SD], in the reply [5T] 
to that comment, and in a recent joint manuscript by the same authors [56j . 
Fluctuations in the mass density arise entirely from fluctuations in the mo- 
mentum and heat fluxes, so that the equations considered in this paper are 
not of a truly hydrodynamic origin. 



4 Macroscopic Fluctuation Theory 

Before turning to macroscopic fluctuation theory, let us summarize some key 
elements of GENERIC. In the GENERIC framework, the building blocks of 
nonequilibrium thermodynamics are formulated to construct time-evolution 
equations for thermodynamic systems. In particular, the hallmark of nonequi- 
librium entropy is to generate irreversible dynamics. For complex fluids, this 
nonequilibrium entropy depends on additional structural variables, such as 
polymer conformation tensors in polymeric liquids [46j . In macroscopic fluc- 
tuation theory [9j (see also tSkZj), the hydrodynamic equations for averages 
obtained as an output from GENERIC are assumed to be given as an input 
to the theory of nonequilibrium thermodynamics because the interest is in 
the fluctuations around the solutions of the hydrodynamic equations. For the 
diffusive system considered in the previous section, it is assumed that one 
knows the continuity equation for the density 

with a mass flux j — j (p) given by the constitutive equation 

j = -D{p)-^ + x{p)-E. (28) 



These two equations correspond to Eq. (|23|) with an additional flux contribu- 
tion resulting from an external field E, where xip) is the mobility tensor. The 
external field can be used to control the evolution of the density profile. The 
hydrodynamic equations (j27p and (|28p are supplemented by boundary condi- 
tions, more precisely, by specifying the chemical potential fi on the boundary 
dA of the domain A. 

In macroscopic fiuctuation theory, the hallmark of nonequilibrium entropy 
is to govern fluctuations. The goal of this approach is to introduce such an 
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entropy entirely in terms of macroscopic concepts and considerations, with a 
validation of the phenomenological ideas by microscopic models Oil]. 

In intuitive terms, the fluctuations around stationary states can be un- 
derstood in the following way (all the details can be found in [6,7,8112]): 
Fluctuations decay according to the macroscopic hydrodynamic equations, 
and they spontaneously emerge according to the time-reversed trajectories 
governed by "adjoint hydrodynamics." Note that the fluctuations are as- 
sumed to satisfy the non-fluctuating boundary conditions and they hence 
depend on the nature of the boundary conditions. Consideration of station- 
ary states is crucial for the analysis of emerging and decaying fluctuations. 
In other words, the construction of adjoint hydrodynamics depends on the 
probability density for fluctuations in the steady state. Only for equilibrium 
states, adjoint hydrodynamics coincides with hydrodynamics; we then have 
time-reversal symmetry. The macroscopic time-reversal behavior of diffusion 
equations has been studied in great detail within the theory of stochastic dif- 
ferential equations in the classical work of Nelson [1^ and, more generally, is 
part of the theory of Markov processes. These time-reversed trajectories have 
been shown to minimize a suitable cost function for fluctuations, which is nat- 
urally identified with the time integral of the extra dissipation rate caused 
by fluctuations. The trajectories of adjoint or time-reversed hydrodynamics 
minimize the work required for a dynamic transition from a stationary initial 
state to a given final state, which can be achieved by suitably chosen external 
forces. As a result of such a variational principle, a nonequilibrium entropy 
or free energy associated with fluctuations can be introduced. 

In terms of equations, adjoint hydrodynamics for the diffusive system is 
governed by the continuity equation p7p with the "reversed" mass flux 



Dip) ■ % Xip) 



^ ^ d 5T{p,p) 
E + 2- ^ ' 



dr Sp 



(29) 



instead of the expression in Eq. ((28|) . The quantity J-{p, p) introduced in 
Eq. ((29l) is referred to as the nonequilibrium free energy of the macroscopic 
state p for a system in the stationary state p [S] . It describes the probability 
distribution for fluctuations Pp{p) around the stationary state p, 

p,(p)_^^ r_^i (30) 



Pp{p) I ^bTq 

The occurrence of a functional T{p,p) in Eqs. (p9)) and pO|) allows for the 
possibility of nonlocal time reversal and fluctuation effects. 

For the actual calculation of JF(p, p), one considers the extra dissipation 
rate as a Lagrangian and passes by Lcgendre transformation to the Hamilto- 
nian formulation. The associated Hamilton- Jacobi equation is then given by 

mm 



Tr-p)-''^P^\Tr-p)--pYr-'^f'\ 



(fr = 0. (31) 



The Hamilton- Jacobi equation (pij) is a convenient starting point for the 
practical calculation of J-'{p,p), be it by veriflcation of guessed solutions or 
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by perturbation theory. It can be rewritten as 



,3 [ 5T r ^ 9 6T ,2 
J OA op or dp 

(32) 

which is particularly convenient when 5J- /Sp vanishes on the boundary so 
that the surface integral on the right-hand side is equal to zero. 

In general, the solution ^(p, p) depends on the dynamic material proper- 
ties xip) E^nd D{p) occurring in the variational problem. In blatant contrast, 
the statistical expressions for the energy E, the entropy S, and hence also for 
the free energy F of GENERIC (as obtained by means of projection opera- 
tor techniques) imply that all these quantities can be expressed in terms of 
a nonequilibrium ensemble and hence do not contain any dynamic material 
information. 

For homogeneous or, in the presence of external fields, inhomogeneous 
equilibrium states p characterized by j{p) — 0, one finds [9] 

Hp, p) - / Ifipir)) - M(r)p(r)] d'r - [ [/(p(r)) - /l(r)p(r)] d'r. (33) 

J A J A 

The occurrence of the Legendre transform of the Helmholtz free energy den- 
sity / is natural because the problem is controlled by the chemical potential 
on the boundaries. Note that, in the first integral in Eq. (1^51) . there actually 
occurs the chemical potential /i(r) associated with the background stationary 
state p{r) in front of p{r). This causes a nontrivially coupled dependence of 
!F on p and p and is important to obtain a functional derivative that vanishes 
on the boundary, 

^ = Mr-)-/i(r), (34) 

which is needed in the verification of the Hamilton- Jacobi equation ()32p . 
together with the Nernst-Einstein relation 



J A ^P 



d_ 
dr 



xip) 



d_5T_ 
dr Sp 



d_ 
dr 



3{P) 



xip)gP=D{p). (35) 

Note that pd'^f/dp^ is the isothermal speed of sound squared and hence 
positive. In terms of the more natural diffusion tensor S introduced via the 
friction matrix (PU)) of the GENERIC approach, the Nernst-Einstein relation 
takes the simpler form 

X{p) = S{p). (36) 

The condition j{p) — for equilibrium is equivalent to the time- reversal 

symmetry jXp) = 3*{p) El- 

The Legendre transform of J-{p^ p) from the density p to its conjugate 
variable, the chemical potential /i, provides the generating functional for 
correlation functions. This approach has been elaborated in Section 4 of [9]. 
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5 Thermodynamics: Averages and Fluctuations 

As we have emphasized in the preceding sections, GENERIC deals with 
the evolution equations for macroscopic averages, whereas macroscopic fluc- 
tuation theory deals with the fluctuations around steady state averages. 
GENERIC is of little interest in formulating hydrodynamic equations be- 
cause the structure of these equations is fixed by the local conservation laws 
for mass, momentum and energy. The more ambitious goal of GENERIC 
is to develop proper and consistent equations for complex fluids with ad- 
ditional slow structural variables, for which there are no conservation laws; 
then, GENERIC provides helpful structural guidance. For complex fluids, 
the equations for the averages need to be formulated first, and a theory of 
fluctuations can only be developed as a second step. Only for conserved quan- 
tities, the balance equations for the averages are so straightforward that one 
can proceed directly to a theory of macroscopic fluctuations. We hence re- 
fer to nonequilibrium thermodynamics of the first and second kind. From a 
statistical mechanics perspective, one could also say that thermodynamics 
of the first kind is associated with the law of large numbers for averages, 
whereas thermodynamics of the second kind is associated with the central 
limit theorem 0. 

There are a number of further associations that should go with the dis- 
tinction of thermodynamics of the first and second kind that go beyond the 
emphasis of a natural sequence of steps. In the case of small Gaussian fluctua- 
tions, the averages of thermodynamics of the first kind are given by the first 
moments, whereas the fluctuations of thermodynamics of the second kind 
are described by the second moments. There is also a correspondence to the 
fluctuation dissipation relations of the first and second kind distinguished by 
Kubo [33]. The fluctuation dissipation relation of the first kind deals with 
the average response of a system to an external perturbation and belongs to 
the field of linear response theory. The fiuctuation dissipation relation of the 
second kind provides information about the second moments of the noise in 
a stochastic description of a system and is deeply linked to the projection 
operator approach. 

At equilibrium, of course, thermodynamics of the first and second kind 
are related by Einstein's theory of fiuctuations (see, for example. Section 10. B 
of [58' or, more generally. Chapter 19 of [13'). Far away from equilibrium, 
an equivalence of thermodynamics of the first and second kind is far from 
obvious. A deeper comparison between GENERIC and the macroscopic fluc- 
tuation theory based on fiuctuation dissipation relations is attempted in the 
subsequent section. Among the other attempts to establish a general theory 
of nonequilibrium steady states mentioned in the introduction, the work of 
Oono and Paniconi (1998) [41] and of Sasa and Tasaki (2006) [62] may be 
classified as thermodynamics of the first kind, whereas the work of Taniguchi 
and Cohen (2007) [51] deals with thermodynamics of the second kind. 

The distinction elaborated in the present section is also useful for the 
discussion of different variational principles. On the one hand, the so-called 
principle of minimal entropy production of linear irreversible thermodynam- 
ics ^57jil4ji32, , that is, within thermodynamics of the first kind, is of limited 
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validity (see, for example, p. 832 of or Section 3.1.5 of On the 

other hand, the minimum dissipation principle assumed in the construction 
of macroscopic fluctuation theory, which was briefly discussed in the para- 
graph before Eq. (P^ . seems to be of more general validity. 

6 Fluctuation Dissipation Relations 

Nonequilibrium thermodynamics of the first and second kind focus on av- 
erages and fluctuations, respectively. In the GENERIC formulation of the 
equations for averages, an important feature is the proper formulation of 
irreversible or dissipative dynamics (after separating it from reversible dy- 
namics). It is hence natural to look at fluctuation dissipation relations in 
order to establish a connection between GENERIC and macroscopic fluctu- 
ation theory. Each of the approaches has offered such relations. 

6.1 From Macroscopic Fluctuation Theory to GENERIC 

In the framework of macroscopic fluctuation theory, there exists a so- 
called nonlinear fluctuation dissipation relation for stationary nonequilibrium 
states. From the equations of SectionlH a version of this nonlinear fluctuation 
dissipation relation is obtained as the sum of Eqs. ([^5]) and (f^^ . 

-,{3+f)^~xip)-^^-^. (37) 

This equation (see also Eq. (2.15) of W) is reminiscent of the irreversible 
contribution to GENERIC in Eq. (fT4|) . however, it is the nonequilibrium 
free energy governing fluctuations that generates the arithmetic mean of 
real hydrodynamic and adjoint hydrodynamic evolution characterizing the 
relaxation and emergence of fluctuations around a given stationary state, 
respectively (for equilibrium states with j = j*, the entropies of GENERIC 
and of macroscopic fluctuation theory hence coincide). The reformulation of 
the nonlinear fluctuation dissipation relation to obtain an equation for the 
relaxation of fluctuations in the top half of p. 645 of [5] is even more similar 
to GENERIC, 

where S is the entropy of macroscopic fluctuation theory and A{p) is a vector 
fleld that conserves the entropy, 

/ ^A{p)d^r^O. (39) 

Equation ([M)) corresponds to the degeneracy ^ of the GENERIC frame- 
work, which is a hallmark of reversibility of the vector field A{p) . GENERIC 
additionally postulates a degenerate Poisson structure to formulate the 
entropy-conserving or reversible vector fleld A{p) in terms of a Hamilto- 
nian. Once more we emphasize that, in spite of striking similarities between 
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Eqs. (fT2|) and ([38|) . the entropy S of the macroscopic fluctuation theory de- 
pends on the underlying steady state and on dynamic material properties, 
whereas the entropy S of GENERIC does not. 

If we neglect the surface term in the Hamilton-Jacobi equation ([32]) and 
rewrite it in the compact form 



Sp 



d^r = 0, (40) 



then we arrive at an alternative interpretation of this equation inspired by 
the GENERIC framework. The free energy T is to be constructed such that 
it generates irreversible dynamics and, at the same time, is conserved under 
the reversible dynamics remaining after subtraction of the irreversible contri- 
bution from the full dynamic equation. This double role of the free energy is 
the reason for the quadratic occurrence of 5T/5p in Eq. (|40|l . Also note that 
the dissipative bracket [J^, T] of Eq. pT|) occurs in the full Hamilton-Jacobi 
equation (|3T|) . 

According to the GENERIC approach, when formed with a constant tem- 
perature To, the free energy F{p) actually generates irreversible dynamics 
and is conserved so that T = F satisfies Eq. ([3D]). However, SF/Sp does not 
vanish on the boundary, as was assumed in order to arrive at Eq. (|40p . The 
difference between the nonequilibrium free energies of GENERIC and the 
macroscopic fluctuation theory hence results exclusively from the boundary 
condition 5F/Sp = 0. The complicated functional dependence of T on the 
underlying steady state p and on dynamic material properties is recognized 
as entirely due to the boundary conditions propagating into our driven dif- 
fusive system to achieve steady state conditions. In short, the fluctuations 
of the steady state system are complicated because no fluctuations of the 
chemical potential on the boundary are allowed. 



6.2 From GENERIC to Macroscopic Fluctuation Theory 

We have seen that the macroscopic fluctuation theory comes with GENERIC- 
type equations, but with a more complicated entropy resulting from bound- 
ary effects. We can now ask the reverse question: Can GENERIC provide 
equations for fluctuations around the solutions of the equations for averages? 
Fortunately, GENERIC is backed up by statistical mechanics and a corre- 
sponding fluctuation dissipation relation. 

The statistical approach to GENERIC is based on projection operator 
techniques and nonequilibrium ensembles (see Section 6.1.2 of [46]), and 
the choice of the ensemble clearly corresponds to an assumption about the 
nonequilibrium fluctuations. This may be considered as a weakness of the sta- 
tistical approach to GENERIC. An analogous situation arises in Einstein's 
theory of equilibrium fluctuations where one needs to argue that the ex- 
tensive quantities fluctuate for given intensive variables, which are flxed by 
a surrounding bath. Landau and Lifshitz are willing to assume an isolated 
system at constant energy to justify a fluctuating temperature '34J. For a 
generalized canonical nonequilibrium ensemble, the probability distribution 
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for fluctuations Pp{p) around a given state p occurring in Eq. (I30|) has ac- 
tually been discussed in a more general context in Exercise 138 of [46]. The 
resulting nonequilibrium free energy associated with the fluctuations of our 
diffusive system can be written as 



T{p, -p) = F[p) - F{p) ~ / /2(r) [p{r) - p{r)] d\ 



(41) 



which is a direct generalization of the equilibrium free energy in Eq. (jSSp . 
Such a simple local equilibrium generalization, and hence the concept of gen- 
eralized canonical nonequilibrium ensembles, would be inappropriate accord- 
ing to the macroscopic fluctuation theory. The deeper reason is that inten- 
sive variables can be prescribed only on the boundaries and not throughout 
the bulk system. Note, however, that macroscopic fluctuation theory in its 
usual form relies on non-fluctuating boundary conditions for the chemical 
potential which may be challenged, at least from a phcnomenological per- 
spective (and one may try to justify them by statistical mechanics). The 
intensive quantities are related to the Lagrange multipliers of the canonical 
ensemble; as such, they are a property of the entire ensemble rather than 
individual fluctuations. Even at equilibrium, the naturally fluctuating vari- 
ables are the extensive quantities, but fluctuations of intensive variables like 
temperature or chemical potential are often introduced through their ther- 
modynamic equations of state (see, for example, Eq. (10.14) of [5S] or § 112 
of [34]; on the other hand, in Chapter 19 of his textbook [13], Callen refrains 
from introducing fluctuations of intensive variables) . 

A more fundamental approach to the fluctuations to be added to 
GENERIC is obtained by using probability densities of the original vari- 
ables as new variables in the projection operator approach (see Section 6.3 
of 06]). It is found that the generalized canonical ensemble on the level of 
probability densities is intimately related to the generalized microcanonical 
ensemble on the level of the original variables. As a result of the projection op- 
erator procedure, one obtains the very natural and appealing Fokker-Planck 
equation. 



dp{x, t) 
dt 



5 
Sx 
, 6 
ox 



r, JE{x) JS{x) 
L{x)^^ + M{x) ^ ' 



Sx 



M{x)—p{x,t) 
ox 



Sx 



p{x,t) 



(42) 



governing the dynamics of fluctuations in arbitrary nonequilibrium systems, 
where all GENERIC building blocks must be evaluated in the fundamental 
microcanonical ensemble. Equation ([42]) consists of a GENERIC drift and 
a superimposed multiplicative white noise |24p43j . The representation of the 
solutions of Fokker-Planck equations in terms of functional integrals and the 
development of fleld-theoretic solution methods can be found in the text- 
book 28J; such a representation is particularly useful for the interpretation 
of Eq. (|42ll in terms of fluctuating trajectories. 



15 



For the evolution of the average of an arbitrary observable A, the Fokker- 
Planck equation (|42|) implies 



f^.(M,£)) + (|A51)+.B(5^M^). (43) 

This equation can be used to evaluate averages and correlations in time- 
dependent and steady state situations. 

We are particularly interested in fluctuations with experimentally observ- 
able consequences, most importantly in scattering experiments [H]. Except 
near critical points, fluctuations are usually very small which, of course, is im- 
portant for the success of thermodynamics. At equilibrium, Einstein's famous 
fluctuation theory hence relies on a second-order expansion of the entropy 
leading to Gaussian fluctuations. Whereas large fluctuations can be treated 
easily and elegantly (see Chapter 19 of [13), there is usually no need to go 
beyond Gaussian fluctuations. In the same spirit, the Fokker-Planck equation 
of GENERIC with fluctuations can handle all problems involving small 
noncquilibrium fluctuations. A special case of the Fokker-Planck equation 
(|42| corresponds to fluctuating hydrodynamics. Whereas the correlations of 
the non-conserved fluxes are short range, in noncquilibrium systems, the 
well-known long range fluctuations of conserved quantities arise even after 
linearization of the noise terms [42] . 

If we are interested in small fluctuations around steady states, we make 
the simplifying assumption that the averages x are characterized by the con- 
dition of vanishing macroscopic time evolution, 

dx ^,_.6E(x) ^^,-.SS(x) „ 

— =L{x)^^+Mix)^^ =0. (44) 

dt ox ox 

In other words, we do not take fluctuation renormalization [251l46j into ac- 
count. It is then natural to assume additive or linearized noise, that is, to 
replace M{x) by M{x) in the last term of Eqs. ^ and From Eq. (03]) 
we then obtain a powerful equation for calculating steady state correlations 
of small fluctuations of x around x, 

{{x„Xj)) + {{x,,xj)) + 2fcBAf„- {x) = 0, (45) 

where ((•,•)) denotes the covariance of two observables and x is the linearized 
right-hand side of the GENERIC evolution equation IT^ . The linearized 
GENERIC evolution involves reversible dynamics and the entropy in addi- 
tion to the friction matrix M . For the diffusion problem, it is the discrepancy 
between S in the friction matrix in Eq. (I20p and D in the evolution equation 
([^5]) that produces interesting correlations. According to Eq. (|35p . this dis- 
crepancy is determined by the behavior of the local equilibrium free energy. 
In hydrodynamics, the reversible convection effects occurring in the first two 
terms of Eq. ([l5]) lead to long range density correlations [32] ■ 

As mentioned at the end of Section [51 the assumption of isothermal con- 
ditions may be subtle and might require a more careful discussion because, 
except in proper limiting cases, a constant Tq is inconsistent with the steady 
state conditions ([44]) and (|45|) . We tried to circumvent this problem by con- 
sidering athermal hard-sphere systems for which a GENERIC formulation 
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without a temperature variable works well. However, a dependable discus- 
sion of the local entropy production and flux is clearly desirable. 



6.3 One-Dimensional Simple Exclusion Process 



In order to illustrate the calculation of long range nonequilibrium correla- 
tions within the GENERIC framework with fluctuations, we consider one- 
dimensional diffusion in the interval [0,^]. We assume a constant diffusion 
coefficient so that the steady state solutions of the diffusion equation ((23|) 
are of the form 



P{r) = Po + {pi - Po) 



(46) 



where po and pi are the values of the density profile p at the boundaries of 
the interval [0, 1]. For the quantity S we assume a particular form of Eq. (|25p . 



1 



P_ 

Pc 



D. 



(47) 



If the parameter pc in this equation goes to infinity, we recover ideal gas be- 
havior. A finite value of pc corresponds to the simple exclusion process of [8] . 
Note, however, that all our calculations are based on the phenomenological 
equation ((44|l and the fluctuation dissipation relation ((45|) . 

If we apply Eq. (|l5|) to Xi = p{r), Xj = p{s) and introduce the covariance 
function 

cov(r,s) = ((p(r),p(s))) , 



then we obtain the equation 



92 



d 



TO + TO cov(r, s) = 2m— p(r) 1 - 



p{r) 



Qj,2 ds^ J ' dr ' \ Pc 

The most general solution of this equation is given by 

p{r) 



^S{r - s). 
or 



(48) 



(49) 



cov(r, s) — mp{r) yi — 
MPI - Pof 



S{r 



PcP 



rs 
T 



(50) 



where (/)(r, s) is a solution of the two-dimensional Laplace equation on the 
square [0, 1] x [0, 1] determined by the desired boundary conditions. We hence 
actually find the previously mentioned freedom of choosing the boundary 
conditions. Nonfluctuating boundary conditions imply 0(r, s) = 0. For r > s, 
we then have the equation 



cov(r, s) ^ ~{pi- Po)- 



m s 
Ipc I 



(51) 



which expresses the negative nonlocal correlations for the simple exclusion 
process given m Eq. (2.3) of [8]. Note that these nonlocal correlations vanish 
in the limit pc 00, that is, for the ideal gas. The exclusion mechanism is 
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crucial for obtaining long range fluctuations in nonequilibrium steady states, 
that is, for pi / po- In the GENERIC approach, there is no need to in- 
troduce a nonlocal entropy functional to obtain the correlations ((50)) form 
purely phenomenological equations. The GENERIC nonequilibrium entropy 
is of the local equilibrium form. The origin of this simplicity lies in the fact 
that the relationship ()35p depends only on local equilibrium thermodynamic 
properties. 

We once more emphasize the possibility of choosing different boundary 
conditions for <j)[r,s). This implies different correlations and hence different 
nonequilibrium free energies according to macroscopic fluctuation theory. 

7 Summary and Discussion 

In the context of an isothermal driven diffusive system, we have compared 
two different approaches to nonequilibrium systems, one based on averages 
and the other one on fluctuations. To emphasize this conceptual difference 
between the approaches, we have classified the GENERIC framework, which 
is focused on the structure of time-evolution equations for averages, and the 
macroscopic fluctuation theory as thermodynamics of the first and second 
kind, respectively. In general, one first needs to find the proper equations for 
the averages before, in a second step, one can discuss the fluctuations around 
them. Deep relationships between the two approaches can be revealed by 
means of the fluctuation dissipation relations that exist in either approach. 

Both approaches introduce nonequilibrium entropies and Helmholtz free 
energies. Whereas the free energy of GENERIC contains only static material 
information, the free energy of the macroscopic fluctuation theory depends 
on transport coefficients. This difference is revealed to be a consequence of 
boundary conditions which, in macroscopic fluctuation theory, are assumed 
to be non-fluctuating. Different boundary conditions lead to different free 
energies. 

The hydrodynamic equation for the decay of fluctuations in the macro- 
scopic fluctuation theory formally possesses the GENERIC structure, and 
the Fokker-Planck equation governing fluctuations within the statistically 
founded GENERIC approach provides the counterpart to the Hamilton- 
Jacobi equation for the nonequilibrium free energy of the macroscopic fluc- 
tuation theory. The steady state moment equation (I45p resulting from this 
Fokker-Planck equation provides a powerful tool for calculating long range 
nonequilibrium correlations. An explicit calculation for a one-dimensional 
model shows that GENERIC can reproduce the long range nonequilib- 
rium correlations of macroscopic fluctuation theory in spite of the different 
nonequilibrium entropies used in the two approaches. The origin of these long 
range correlations for a model with dissipative mass flux, which is not truly of 
hydrodynamic origin, lies in the local equilibrium thermodynamic properties 
and their position dependence resulting from a nonuniform density proflle. 

It is well known that the different nonequilibrium ensembles used in the 
projection operator derivation of GENERIC may at best be equivalent for the 
resulting average equations but certainly not for the description of fluctua- 
tions. Macroscopic fluctuation theory suggests that the generalized canonical 
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ensemble does not represent fluctuations in a boundary-driven system in an 
appropriate way. However, this conclusion relies on the debatable assump- 
tion of non-fluctuating boundary conditions. A more systematic projection 
operator based theory of fluctuations, even for time-dependent nonequilib- 
rium systems, is given by the Fokker-Planck equation (j42|) . A detailed dis- 
cussion shows the possibility of choosing from a variety of different bound- 
ary conditions. An even more illuminating and conclusive comparison be- 
tween GENERIC and macroscopic fluctuation theory might be obtained for 
flow-driven systems, where the driving mechanism is felt via a velocity field 
throughout a bulk system rather than at the boundaries only and a reliable 
complete set of hydrodynamic equations is used. 

This paper is a plea to look for simplicity in beauty and generality rather 
than in propitious special cases. Eventually, fully consistent results should 
be achieved. 
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